arXiv:1509.04787vl [cond-mat.mtrl-sci] 16 Sep 2015 


Substrate effects and diffusion dominated roughening in CU 2 O electrodeposition 


1. S. Brandt^, V. C. Zoldan^, V. Stenger^, C. C. Pla Cid^, A. 

A. Pasai’(“), T. J. Oliveira^’^'') and F. D. A. Aarao Reis^’^^) 

^ Departamento de Fisica, 

Universidade Federal de Santa Catarina, 

88040-900, -^76, Flonanopolis, SC, Brazil 
^ Departamento de Fisica, 

Universidade Federal de Vigosa, 

36570-000, Vigosa, MC, Brazil 
® Instituto de Fisica, Universidade Federal Fluminense, 

Avenida Litoranea s/n, 24210-340 Niteroi RJ, Brazil 

(Dated: September 17, 2015) 

Cuprous oxide (CU 2 O) films from 25 nm to 1500 nm were electrodeposited on n-Si(lOO) and Ni/n- 
Si(lOO) substrates from aqueous solution at room temperature. X-ray diflfraction and transmission 
electron microscopy imaging show that the CU 2 O structure and morphology is strongly aflfected by 
the substrate choice, with V shape and U shape columnar growth on n-Si(lOO) and Ni/n-Si(100), 
respectively. Atomic force microscopy reveals the presence of rounded grains at the surface in both 
cases. Anomalous and normal roughening are observed in films grown on n-Si and Ni, respectively, 
but estimates of scaling exponents are not conclusive. On the other hand, the distributions of 
local heights, roughness, and extremal heights show good agreement with those of the fourth order 
linear stochastic equation of Mullins and Herring (MH). Thus, surface dynamics in both systems is 
dominated by diffusion of adsorbed molecules, with no large scale effect of possible inhomogeneities 
in mass flux from the solution or in reaction and adsorption rates. In growth on n-Si substrates, the 
noise amplitude of the MH equation increases in time as while the coefficient of the curvature- 
related term is time-independent. Step edge energy barriers restrict the mass flux across grain 
boundaries, thus a broad size distribution of initial grains leads to coarsening of the larger ones. 
This explains their V shape in the thickest films and establishes a connection with the anomalous 
roughening. These effects are reduced in films grown on Ni/n-Si, which initially have much larger 
grains with narrower size distributions and, consequently, smaller fluctuations in coarse grained 
growth rates. Thus, despite the relevance of electrochemical conditions for CU 2 O films to grow and 
their influence on crystallographic orientation, large scale surface features are determined by physical 
properties of the material and its interactions with the substrate, with a universal microscopic 
dynamics similar to vapor deposition. 


I. INTRODUCTION 

CU 2 O (cuprous oxide) is a p-type semiconductor that 
has recently attracted the attention of experimental and 
theoretical groups due to its potential for application in 
metal base transistors [T], spintronic [2], photocathode 
for water splitting [3] , electrochemical supercapacitors [4] 
and catalysts [5], and for light harvesting [6]. The de¬ 
position of CU 2 O layers has been achieved by different 
techniques, such as pulsed laser deposition [7], magnetron 
sputtering [8] , copper oxidation [9] , radical oxidation [10] , 
and electrodeposition m- Electrodeposition is a very 
versatile technique, allowing to modify many structural, 
optical, and electrical properties by the control of the 
deposition parameters [IMl]. Moreover, electrodepo¬ 
sition appears as a cost effective method to the prepa¬ 
ration of metal and semiconductor thin films for device 
applications. When compared with physical deposition 
processes, such as sputtering, electrodeposition has the 
particular feature of diffusive mass transport of species 
from the electrolyte to the growing surface [TO] . 

A comprehensive understanding of the electrochemi¬ 
cal growth mechanisms, taking into account the mass 


transport at the electrode surface, plays a vital role to 
the production of films with the desired properties. One 
example is the technological requirement for mesoscopic 
layers to be used in photocathode applications, in which 
the grain shape and the exposed crystalline faces need to 
be controlled to improve stability and efficiency [IT]. On 
the other hand, the study of surface topography of de¬ 
posited layers helps to discriminate the basic mechanisms 
of growth dynamics from the particular features of each 
set of physico-chemical deposition conditions [THHUj- A 
crucial condition may be the substrate, which affects the 
early stages of island formation and growth and, possibly, 
the subsequent film dynamics. 

Bearing in mind the above mentioned arguments, the 
present work aims to investigate the structure of electro- 
chemically grown CU 2 O films on top of a semiconducting 
and monocrystalline substrate and on top of a metal¬ 
lic and polycrystalline substrate, viz. n-Si(lOO) and Ni 
evaporated on n-Si(lOO). Recently, our group showed 
that these two substrates are suitable for morphological, 
structural and optical studies on CU 2 O [13]. Moreover, 
the Cu 20 /Ni structure may be employed on photochem¬ 
ical [53] and spin transport [23] investigations. 
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Raman spectroscopy, X-ray diffraction (XRD), and 
transmission electron microscopy (TEM) measurements 
are carried out to characterize film composition and 
structure. Atomic force microscopy (AFM) images are 
used to study kinetic roughening of the film surface, thus 
helping to disclose the main mechanisms of the growth 
dynamics. Films directly grown on the semiconductor 
substrate have a V shape columnar structure, preferen¬ 
tial growth in [ 100 ] direction, and show anomalous rough¬ 
ening, while those grown on the Ni layer present an U 
shape columnar structure, change the dominant crystal¬ 
lographic orientation during the growth, and have nor¬ 
mal roughening. A deeper analysis of surface morphol¬ 
ogy shows that the dominant mechanism in the CU 2 O 
growth is surface diffusion of adsorbed molecules, with 
uniform incoming flux. Step edge energy barriers explain 
the anomaly of the former films, connecting this feature 
to the grain shape. The universal coarse-grained growth 
dynamics is similar to what is observed in many vapor 
deposited films, despite the relevance of electrochemical 
conditions to local reaction and adsorption processes. On 
the other hand, the interaction between the film and the 
substrate, which is frequently restricted to island nucle- 
ation and growth, here is shown to have drastic conse¬ 
quences to the structure of the thickest CU 2 O films. 

The rest of this work is organized as follows. In Sec. 
im the experimental procedure for growth and character¬ 
ization of the CU 2 O films is described. In Sec. [ml the 
experimental results are presented and analyzed in the 
light of kinetic roughening concepts. Section |TV| presents 
a discussion of the results and our conclusions. 


II. EXPERIMENTAL 

CU 2 O films were deposited on ( 100 ) monocrystalline 
n-type silicon (resistivities in the range of 6 - 9 fl-cm) 
without and with a cap layer of 50 nm of evaporated 
Ni. Before CU 2 O electrodeposition or Ni deposition, the 
silicon substrates were immersed in HF 5% for 20s to 
remove silicon native oxide from the surface. The rough¬ 
nesses of the n-Si substrate and of the Ni layer surface are 
respectively 0.2 nm and 1.3 nm. Ni thin films deposited 
on Si(lOO) had preferential growth in the [111] direction, 
which was checked by XRD measurements. 

The electrochemical experiments were conducted in a 
conventional three electrode cell connected to a potentio- 
stat Autolab PGSTAT30. The counter electrode was a Pt 
foil and the reference electrode was a Saturated Calomel 
Electrode (SCE). The electrolyte, prepared from analyt¬ 
ical grade reagents and deionized water (resistivity of 18 
M fi-cm), contains 0.4 M CUSO 4 and 3.0 M lactic acid, 
with the pH adjusted to 10.00 by adding a solution of 5.0 
M NaOH [ni- 

The deposits were obtained via potentiostatic experi¬ 
ments. Samples were deposited at 25 °C for a deposition 
potential of —0.5 V vs SCE, where the only varied pa¬ 
rameter was the thickness H of the samples for the 2 


different substrates, n-Si(lOO) and Ni/n-Si(100). 

The efficiency of the CU 2 O electrodeposition process 
was checked by Rutherford backscattering in a previous 
work and values of 90% were obtained The sam¬ 
ples were characterized by various techniques including 
Raman spectroscopy (inVia, Renishaw), XRD (XPERT, 
Philips), TEM (JEM-2100, JEOL) and AFM (Pico-SPM, 
Molecular Imaging Corporation). Raman spectra were 
obtained from 514.5 nm wavelength excitation (argon 
ion laser). The AFM images were acquired in contact 
mode in air at room temperature and the values of the 
root mean-square deviation of the local surface height, 
W{l,H), are representative of at least three images of 
different regions of the surface. Three different regions 
for sample were measured with scan sizes of 1 x 1 , 2 x 2 
and 5 x h/rm for each region and number of pixels per 
line ranging from 256 to 1024, respectively. From these 
measurements we have concluded that for 2 x 2fim image 
size the W{l,H) reaches the saturation value. In addi¬ 
tion, before starting measurements on CU 2 O samples, the 
performance of the AFM tip was checked by imaging Au 
thin films with known surface. 


III. RESULTS 

A. Characterization of film composition and 
structure 

The copper oxide phase formed in the electrodeposited 
thin films was analyzed by Raman spectroscopy and 
XRD. The measurements were done on layers with thick¬ 
ness between 150 and 1500 nm deposited on n-Si(lOO) 
and Ni/n-Si(100). 

Raman spectroscopy results of films with thickness 
1250 nm are shown in Fig. The predominance of 
CU 2 O phase is evident. Contributions of CU 4 O 3 and pos¬ 
sibly CuO phases are also observed. The peaks in Fig. 

were addressed to the respective phase and Raman 
vibration mode based on Ref. [26]. Despite the obser¬ 
vation of CU 4 O 3 and CuO Raman peaks, these phases 
are not present in XRD Bragg-Brentano patterns of the 
same films, shown in Fig. [^, probably due to the higher 
surface sensitivity of Raman measurements. CU 4 O 3 and 
CuO phases are likely to be formed by oxidation of CU 2 O 
top layers. Since the layers are mostly of CU 2 O, this work 
will be concerned to the description of the electrochemi¬ 
cal growth of this phase. 

Fig. ^ also shows the XRD patterns of CU 2 O films 
with 1250 nm of thickness deposited on Ni/n-Si(100). 
The peaks are also in the expected positions for the CU 2 O 
crystal. Cu 2 O/n-Si( 100 ) samples display strong prefer¬ 
ential growth in [ 100 ] direction, following the orientation 
of the substrate. On the other hand, CU 2 O film grown 
on top of Ni/n-Si(100) substrate is composed by grains 
oriented in [100], [110], [111], and [311] directions. 

The evolution of CU 2 O growth on Ni/n-Si(100) is ex¬ 
plained by Ref. [T3| as follows: (i) the CU 2 O films 
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FIG. 1. (Color online) a) Raman spectroscopy and b) XRD 
results for CU 2 O films of 1250 nm electrodeposited on n-Si 
and Ni substrates. For each Raman peak are indicated the 
possible related crystalline phases and, between parentheses, 
the Raman vibrational mode. XRD characterization was car¬ 
ried out in Bragg-Brentano configuration. The inset in b) 
shows the CU 2 O portion growing in [100] and [111] directions 
as a function of the film thickness. 


show initial [111] growth due to a better coupling 
of the Cu20(lll)/Ni(lll) interface compared to the 
Cu2O(100)/Ni(lll) one; {ii) at electrolyte pH of 10.00 
the [100] growth is favored, the initial [111] growth starts 
to be suppressed by the [100] one and at a specific thick¬ 
ness a major part of the CU 2 O deposit will be composed 
by [100] grains. 

These growth steps are confirmed by the inset of Fig. 
[^, which shows the fraction of the material growing in 
[100] and [111] directions as a function of film thick¬ 
ness for both substrates. While CU 2 O on n-Si(lOO) is, 
independently of film thickness, oriented in [100] direc¬ 
tion, CU 2 O on Ni/n-Si(100) grows up to 750 nm with 
~ 60% of its grains oriented in [111] direction, but only 
~ 30% remain with this orientation in larger thicknesses. 
The initial growth of CU 2 O on Ni(lll) is thermodynami¬ 
cally controlled by the Cu20(lll)/Ni(lll) coupling, like¬ 



FIG. 2. (Color online) Cross section TEM images of CU 2 O 
thin films electrodeposited respectively on a) n-Si and b) 
Ni/n-Si. Different layers in each image are labeled as fol¬ 
lows: 1 CU 2 O film, 2 n-Si substrate, 3 Ni film (substrate) and 
4 the glue used on sample preparation. In a) and b) discon¬ 
tinuous red lines highlight the V shape and U shape of the 
CU 2 O columns, respectively. 


wise as previously reported for Cu20(lll)/Au(lll) and 
Cu20(100)/Au(100) interfaces [12]. However, the crys¬ 
tallographic orientation transition at 750 nm is related 
to the oxygen concentration in electrolyte at pH 10.00, 
which favors the [100] growth [13] . 

The structure of CU 2 O films was checked by cross sec¬ 
tion TEM images. Figure shows an image obtained 
from electrodeposited CU 2 O film on n-Si(lOO). A high 
density columnar microstructure without voids is ob¬ 
served. These columns have a V shape as stressed by 
discontinuous red lines, exhibiting a grain width increase 
as a function of film thickness. On the other hand, as seen 
in Fig. CU 2 O films deposited on the Ni layer under 
same conditions (solution, temperature and deposition 
potential) displayed a columnar microstructure with U 
shape, as indicated in Fig. and grain width roughly 
constant as a function of growth time. 

Figure shows 2 x 2/xm^ AFM images of the films 
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FIG. 3. (Color online) 2/im x 2fim AFM images of the surface 
morphology of CU 2 O layers electrodeposited on: a), b) n-Si; 
and c), d) Ni/n-Si. Upper row is for thickness of 250 nm, 
while lower row is for 1500 nm. Bright areas are elevations 
and dark ones are deep regions. 

with thicknesses of 250 and 1500 nm. A common fea¬ 
ture in both substrates is that the film surface presents 
a granular aspect. In films deposited on n-Si, there is a 
remarkable increase of the grain size from the thinnest to 
the thickest film. The columnar structure with V shape 
shown in TEM images and this increase in grain size sug¬ 
gest a coarsening process, with part of the initial columns 
growing and enlarging up to 1500 nm at the expenses of 
neighboring ones. On the other hand, a mild increase of 
grain size is observed in the films grown on Ni/n-Si, which 
is in agreement with TEM images showing columns in U 
shape. 


B. Dynamic scaling of surface roughness 

The surface roughness W is the most used quantity to 
characterize height fluctuations. For a given surface un¬ 
der a process of kinetic roughening, W depends on the 
lengthscale in which fluctuations are measured. More¬ 
over, given a scanning box with lateral size I, the fluctu¬ 
ation in film height inside the box also depends on the 
box position. 

For these reasons, the box roughness is defined as the 
root-mean-square (rms) fluctuation of the film height (h) 

inside a box at a given position, i.e, w = \J— h . The 
average local roughness W{I, H) = {w) of a film of thick¬ 
ness H in the lengthscale I is the average of the box 
roughness w among all possible box positions. W{l,H) 



FIG. 4. (Golor online) Roughness W{l,H) as a function of 
the box size I measured from AFM images of a series of elec¬ 
trodeposited GU 2 O films for thickness H from 25 to 1500 nm. 
In a) and b) the substrate is n-Si and Ni/n-Si, respectively. 

is hereafter simply called roughness. 

Using the heights from AFM images, W{l,H) is cal¬ 
culated for CU 2 O layers grown on both substrates. The 
results are shown in Figs. and[^ for a range of thick¬ 
nesses H varying from 25 to 1500 nm. The structural 
changes induced by substrates n-Si and Ni/n-Si are re¬ 
flected in topographic differences in the growing surface 
of the films, confirming the visual inspection of AFM im¬ 
ages. 

The roughness W{I, H) of films grown on Ni/n-Si show 
normal dynamic scaling because it is approximately time- 
independent for small box size I and saturates at val¬ 
ues increasing with the thickness H. Deviations for the 
thinner samples (25 to 75 nm) are characteristic of the 
initial island growth and coalescence, indicating that an 
extended film is formed between 75 nm and 150 nm of 
thickness. For this reason, only results for H ^ 150 nm 
will be analyzed below. 

On the other hand, the oxide layers grown on n-Si have 
anomalous scaling, since the roughness in small boxes 
shift to larger values with increasing thicknesses [27U30] . 
The presence of a single dominant orientation during the 
growth (shown by XRD) and the coarsening of grow¬ 
ing columns from the smallest to the largest thicknesses 
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FIG. 5. (Color online) a) Local roughness for box size I — 40 nm (W*) and b) global roughness (144) as functions of the film 
thickness H. Rescaled correlation functions r(/)/r(0) against the length I for c) n-Si and e) Ni/n-Si substrates. Their first 
minima (the correlation length ^) are shown in d) and f) as functions of H. Continuous lines are power law fits of the data and 
dashed lines have the slopes indicated. 


(suggested by TEM and AFM images) indicate that the 
scaling anomaly is an intrinsic feature of the surface dy¬ 
namics in this system. 

The data presented in Fig. |^for the thickest films with 
I < 100 nm were fit to the relation W ^ l°‘‘. It gives es¬ 
timates of the local roughness exponents ai = 0.90 and 
0.91 in n-Si and Ni/n-Si substrates, respectively. These 
values are consistent with the interpretation of granular 
effects in Refs. [31] and [32|, which predict an exponent 
near 1 if the maximal box size is smaller than the grain 
size. Inspection of the AFM images (Fig. supports 
this hypothesis. These estimates of ai must not be in¬ 
terpreted as true exponents of dynamic scaling, but as 
exponents representing the grain shapes |32] . 

The local roughness W* = W{1*, H) is measured for a 


fixed small box size I* = 40 nm. It is shown in Fig. as 
a function of the thickness H. For films deposited on Ni, 
W* fluctuates around a constant value as the thickness 
increases, from small values to the largest ones. This is 
consistent with normal scaling. On the other hand, for 
films deposited on n-Si, W* shows a power law increase 
from small thicknesses {H = 150 nm) to the largest ones 
{H = 1500 nm). This is the main evidence of anomalous 
scaling in these hlms. The fit of the data in Fig. ii 
by the scaling relation W* ~ H - m gives an anomaly 
exponent k = 0.49 (k is also called local slope exponent 
in theoretical works and local growth exponent jSioc in 
experimental works EilSOj). 

The saturation values of the roughness in Fig. give 
the global roughness 444, which characterize fluctuations 
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of the whole surface at a given thickness. In Fig. [ 1 >, hF. 
is shown as a function of H for both films. The fits to the 
scaling relation ~ for H ^ 150 give the growth 
exponents /3 = 0.84 and (3 = 0.33 for films deposited on 
n-Si and Ni/n-Si substrates, respectively. 

Since the growth on Ni/n-Si substrates has normal 
scaling, we search for a model of local stochastic growth 
equation, which may disclose the main mechanisms of 
the surface dynamics [m HD]. The above estimate of 
/3 is larger than the values provided by the equations 
of second and fourth orders, which are in the range 
[0,0.25] [TD|. However, Fig. shows that a slope close 
to /3 = 0.25 cannot be discarded for large thicknesses. 
This is close to the exponents of the Kardar-Parisi-Zhang 
(KPZ) class [35], /3 « 0.24 |331|3S], and of the Mullins- 
Herring (MH) class [55], P = 0.25. 

Information on surface correlations can be ob¬ 
tained from the spatial correlation function r(/,t) = 

h{x + l,t)h{x,t)'^ at time t, with h = h — {h). The 

rescaled functions T {I, t) of the oxide films respectively 
grown on n-Si and Ni/n-Si substrates are shown in Figs, 
i: an as a function of the distance I, for several thick¬ 
nesses. A reliable estimate of the correlation length ^ is 
the minima of F (Z, t) [37] [35], which is shown in Figs, 
and [If as a function of the thickness H. 

The estimates of ^ are close to those obtained by inter¬ 
section of straight line fits of the roughness in the growth 
and saturation regimes, which is a usual method in the 
analysis of experimental data. Ref. [39] shows that this 
method leads to estimates of ^ of the same order of the 
average grain size estimated by inspection of microscopy 
images [33]. However, an accurate calculation of the aver¬ 
age grain size is more complicated [38j and not performed 
here. 

The scaling ^ ^ ~ is expected for large thick¬ 

nesses [MHO]. Linear fits of the data for H ^ 500 nm 
give estimates Ijz = 0.32 and 0.21 (dynamic exponents 
z = 3.1 and 4.7) for growth on n-Si and Ni/n-Si sub¬ 
strates, respectively. These exponents are very differ¬ 
ent from the KPZ value z « 1.6 [DU [40]. On the other 
hand, the Figs. and [^ show that the slopes of the 
log^ X logiL plots for the largest thicknesses are close 
to the MH value 1/z = 0.25 (the corresponding slope of 
KPZ scaling 1/z = 0.625 is also shown for comparison). 

The MH equation is a model for growth dominated by 
surface diffusion of adsorbed species [55] : 

8h 

— = -KV‘^h + r]{x,t), (1) 

where h (a;, t) is the interface height at substrate position 
X and time t, AT is a constant, and ry is a Gaussian (non¬ 
conservative) noise with covariance amplitude D [191136] . 
A constant external flux was omitted from Eq. Q, which 
corresponds to a spatially uniform and time-independent 
adsorption rate. Thus, the local growth rate of the MH 
model is not affected by diffusion in solution, for instance 
due to the presence of diffusive layers, shadowing effects 
etc. 


On the other hand, the KPZ equation is a model for 
interface growth dominated by surface tension and by a 
nonlinear effect of the local slope [33] : 

^ +r]{x,t), (2) 

where v and A are constants. Distinguishing these pos¬ 
sibilities is an important step to understand the mi¬ 
croscopic growth dynamics of the oxide films. The 
KPZ model also assumes spatially uniform and time- 
independent adsorption rate, with no effect of diffusion 
in solution to the local growth rate. 

C. Scaling of distributions 

We now turn to the analysis of distributions of heights 
and of local quantities, which was recently shown to be 
a more powerful tool to the study of kinetic roughening 
of thin films [371 [m SI]. 

The height distribution P{h) is defined as the proba¬ 
bility density of the height h, so that P{h)dh gives the 
probability of finding a height in the range [h,h dh]. 
Figures]^ and|^, and Ed and w show the scaled height 
distributions of the oxide films with large thicknesses, re¬ 
spectively grown on n-Si and Ni/n-Si substrates. They 
are compared with distributions of the KPZ and MH 
classes in the growth regimes. The KPZ distribution was 
numerically calculated in Refs. [43] and [44]. For the MH 
class, the distribution is Gaussian, which was confirmed 
by direct integration of the MH equation [Eq. [I] an d by 
simulation of the large curvature model (LGM j[3S]. 

The experimental data in Figs. [eji, [6 |d, [^ and w 
are closer to the Gaussian/MH curve. The symmetric 
shape clearly contrasts with the asymmetry of the KPZ 
curve. Indeed, the skewness S and kurtosis Q of the 
P{h) distributions, shown in Figs. and |^, drops to 
zero as the thickness increases (consistently with Gaus¬ 
sian/MH), while the asymptotic KPZ values are S « 0.43 
and Q « 0.34 [431|44]. 

Another important conclusion emerges from Figs. [§i 
an d[^: since the scaled height distributions are the same 
for films deposited on n-Si and on Ni/n-Si, the main 
physico-chemical processes responsible for their rough¬ 
ening are the same. 

The distribution P{w 2 ) of the squared box roughness 
W 2 = is defined so that P{w 2 )dw 2 is the probability 
that the squared roughness in a given box is in the range 
[w 2 ,W 2 dw 2 ] [IDIIIT]- Figures and[^ compares the 
scaled roughness distributions of the oxide films grown on 
both substrates and that of the MH class in the growth 
regime (< 7^,2 is the rms fluctuation of UI 2 ). The size of the 
box in which W 2 is measured ranges from 100 nm to 300 
nm. The MH curve is calculated by integration of the 
MH equation in the growth regime. 

In Figs. and [^, the collapse of the scaled distri¬ 
butions for films grown in different substrates is clear, 
which is additional support to the proposal that their 
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(h - <h>)/W 




FIG. 6. (Color online) Scaled height distributions for a) and c) n-Si, and b) and d) Ni/n-Si substrates. In a) and b) the data 
are presented in semi-log scale and in c) and d) in linear scale. The skewness S and kurtosis Q of the height distributions as 
functions of the thickness H are shown in e) and f) for n-Si and Ni/n-Si, respectively. 


roughening is governed by the same physico-chemical 
processes (despite the very different structures of those 
films). Good agreement with the MH curve reinforces the 
proposal that diffusion of adsorbed species is the main 
mechanism of the surface dynamics. 

We also analyzed the distribution P{m) of the maxi¬ 
mal relative height m = — (h) in the growth regime. 

Here, hm and {h) are, respectively, the maximal and the 
average heights measured inside each box position. The 
distributions for both films are shown in Figs. & and 
showing good data collapse among them and with 
the MH curve. 

Following the procedure described in Ref. [48] , we also 
calculated the local surface slope 6 using the AFM im¬ 
ages. The local slope images corresponding to the height 
images shown in Fig. are displayed in Fig. [^ They 


clearly show that larger values for 9 are obtained in in¬ 
tergranular regions, while small 9 are related to flat areas 
near the top of CU 2 O grains. 

The distributions of local slopes are shown in Figs. [^ 
and for Cu 20 /n-Si and Cu 20 /Ni films, respectively. 
The position 9 of the peaks of those distributions are 
displayed in Fig. as a function of thickness. For films 
grown on n-Si, the distribution significantly changes as 
the thickness increases, with an increase in 9. This is 
typical of anomalous scaling |27l |28| and is qualitatively 
consistent with the local roughness increase shown in Fig. 

On the other hand, for films grown on Ni/n-Si, av¬ 
erage local slopes decrease with the thickness, saturating 
in the thicker films. This is consistent with normal scal¬ 
ing. Comparison of local slope distributions with those 
of theoretical models is not possible because there are 












FIG. 7. (Color online) Scaled a) squared roughness and b) maximal relative height distributions for both substrates and several 
thickness H and box sizes I (measured in pixels of 256 x 256 images). The same data are shown in linear scale in c) and d). 


significant thickness effects and the present concept of lo¬ 
cal slope differs from the definitions in integrated growth 
equations or in lattice models. 

The MH scaling was formerly observed in electrochem- 
ically deposited Prussian Blue films [JS]. In that case, the 
anomalous roughening was generated by time increasing 
adsorption rates. 


D. Extended analysis of MH scaling 

The exact solution of the MH equation [20l [50] allows 
an extension of the interpretation of the experimental 
data, particularly those of films grown on n-Si substrates. 

The MH roughness W{l,H) is expected to scale sim¬ 
ilarly to the height-height correlation function; in terms 
of box size I and time t we have ^D] 

(3) 

where T is a scaling function [F(x) « const for small x] 
and 

^={2Ktf/\ (4) 

Recall that H is proportional to the time t. 

The long time (large thickness) scaling of the cor¬ 
relation length gives exponent z near the MH value 
(1/z ~ 1/4). Thus, Eq. [^indicates that K is approx¬ 
imately time-independent. The coefficient K in Eq. 


represents the relation between surface diffusion coeffi¬ 
cients and local surface geometry of a given material. 
This relation is not expected to change in the course of 
the (constant temperature) deposition. 

For I «C C; Eq. suggests that the loc al rou ghness 
W* (measured for fixed I* =^40 nm; Sec. HI B) scales 
as W*^ ~ ^IniJ. Figure shows W*"^] IniJ versus 
H for both films. The approximately constant value of 
IT*Ini? for the films grown on Ni/n-Si indicates that 
D is also time-independent in that case. However, the 
scaling in films grown on n-Si and the time-independence 
of K gives D ~ Since D is the amplitude of time and 
spatial fluctuations of deposition rate, this means that 
those fluctuations increase in time in the films grown on 
n-Si. This leads to the anomalous scaling in those films. 

The above increase of D is sufficiently slow, thus the 
original solution of the MH equation can be consis¬ 
tently used with the time-dependent form of that pa¬ 
rameter nnnsni- This feature was already illustrated in 
other continuous and lattice models with time dependent 
couplings [nil Eg showing anomalous scaling. 


IV. DISCUSSION AND CONCLUSION 

The scaling of distributions of local heights, maximal 
heights, and roughness provide strong evidence that the 
roughening in electrodeposited CU 2 O films is governed by 
the MH equation Q. The noise amplitude increases in 
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FIG. 8. (Color online) Slope images of CU 2 O films on n-Si, a) and b), and Ni/n-Si, c) and d), evaluated from height images 
shown in Fig. In the first and second row CU 2 O films have 250 and 1500 nm of thickness, respectively. The color bar 
indicates the surface slope 8, which can vary from 0 to 90°. Slope histograms for e) Cu20/n-Si and f) Cu20/Ni structures, and 
g) mean surface slope as a function of the film thickness. 



FIG. 9. (Color online) W*’^/ In H as a function of the thick¬ 
ness H for films grown on n-Si (triangles) and on Ni/n-Si 
(squares). The solid line, with slope 0.8, is a least squares fit 
of the large thickness data for films grown on n-Si. 


time in the equation associated to the films grown on n-Si 
substrates, corresponding to the scaling anomaly. This 
analysis clearly advance over the calculation of scaling 
exponents, which however are also compatible with the 
MH ones. 

The MH equation was proposed for roughening in 
molecular beam epitaxy m, in which a crystalline film 
slowly grows by random adsorption of atoms/molecules 


from vapor and surface adatoms diffuse and preferentially 
bind to high coordination (low energy) sites. Thus, in the 
electrochemical deposition of CU 2 O films, those results 
indicate that this interplay of surface diffusion and de¬ 
position noise is also the main ingredient for roughening 
under the present growth conditions. 

The MH equation considers that adsorption rates are 
spatially uniform and do not depend on time. This sug¬ 
gests that electrochemical conditions do not determine 
the main dynamic mechanisms for roughening, although 
they probably affect the values of the parameters in the 
associated MH equation. Since the same electrochemical 
conditions were used for growth of films on both sub¬ 
strates, we are led to the conclusion that their different 
structures are consequence of the substrate-film interac¬ 
tion. The different interactions in different substrates 
determine the island nucleation and growth and has con¬ 
sequences in the thickest films, despite the universal MH 
dynamics 

The connections between the microscopic dynamics 
and the structure of the growing films can be derived 
from our results, as follows. 

In systems with diffusion-dominated dynamics, 
molecules at a surface grain randomly move, with rates 
depending on local energy barriers. At first approxima¬ 
tion, those barriers depend on coordination numbers and 
increase at step edges. For an adatom or ad-molecule to 
go down a step, the additional activation energy, known 
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as Ehrlich-Schwoebel (ES) barrier [^[^, reduces the 
downward flux near the step edge, favoring nucleation of 
new atomic terraces over the ones previously formed m- 
The effect of ES barriers is enhanced at grain boundaries, 
constraining the mass flux between the grains if height 
differences are large. These features are observed in films 
of several materials grown by vapor methods. The above 
results strongly suggest to extend this interpretation to 
CU 2 O electrodeposition. 

Figure shows that films grown on n-Si substrates 
have an initially small average grain size, but with a 
broad size distribution. This is confirmed quantitatively 
by the small values of the correlation length in Fig. w 

« 80 nm in 75 nm thick films). Thus, the islands nu¬ 
cleating at that substrate also had large size fluctuations. 
This may be a consequence of large spatial fluctuations of 
the adsorption rate during the first stages of the growth. 

When an extended film is formed, the local adsorp¬ 
tion rates may become uniform, particularly if the dom¬ 
inant crystallographic orientation does not change, as 
shown in Sec. |III A[ see also inset of Fig. [^. How¬ 
ever, larger grains have larger numbers of terraces and 
steps, consequently the mass flux from their tops to their 
boundaries is small. With restricted mass exchange in 
grain boundaries, the large grains increase in width and 
cover the smaller ones. From a coarse grained point of 
view, the fluctuations in the adsorption rate continue to 
increase, which qualitatively explains the time increase 
of the parameter D in the corresponding hydrodynamic 
model (the MH equation). 

Figure]^ shows that the shape of the slope distribu¬ 
tion is the same while its average value moves to the right 
during the film growth. This means that local slopes in¬ 
crease inside the grains and at their boundaries, consis¬ 
tently with the above interpretation. 

On the other hand. Fig. n shows that films grown 
on Ni/n-Si substrates have larger initial grains with a 
much narrower size distribution (see also insets of Figs. 

and E)- This is related to small fluctuations in local 
adsorption rates, a feature that is maintained during the 
growth and leads to normal roughening. 

AFM images also reveal that the grain top surfaces are 
similar in the films grown in both substrates, thus the lo¬ 
cal (microscopic) adsorption rates at those surfaces are 


expected to be the same. In growth on n-Si substrates, 
the difference in the net growth rate of individual grains 
have to be explained by the restricted mass flux across 
grain boundaries, not by intra-grain features. If some 
grains with large slopes are initially formed, they may 
grow faster at the expenses of the neighboring ones. This 
establishes a connection between the V-shaped grains 
and the anomalous roughening. 

Unusual scaling features are frequently found in 
diffusion-dominated growth in the presence of step energy 
barriers [H]. For instance, a small barrier in Fe/Fe(100) 
epitaxy leads to mound formation at the film surface, 
with skewed height distribution and f3 = 0.2 m, which 
are characteristics of the non-linear molecular beam epi¬ 
taxy equation Moreover, mound steepening in 

Ag/Ag(100) epitaxy at room temperature occurs with an 
effective exponent /3 much larger than 0.5 in thicknesses 
ranging from 100 to 1000 monolayers [58] (/3 = 0.5 is a 
feature of completely uncorrelated growth (HI). 

We conclude that, although control of electrochemical 
conditions is essential to enable growth of CU 2 O films and 
to determine the crystallographic orientation, our results 
show that large scale surface features are determined by 
physical properties of the material and the substrate, 
with consequences on the internal film structure. Rough¬ 
ening is governed by the interplay of deposition and ad- 
species surface diffusion, similarly to vapor deposition, 
with no significant influence of inhomogeneous mass flux 
from the solution (if it exists). However, substrate-film 
interactions have a strong effect on island nucleation and 
growth, thus the structure of the tickest films is dras¬ 
tically affected by these initial conditions, despite the 
universal (MH) roughening dynamics. 
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